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Abstract 



We present a novel model to calculate vertical transport properties such as 
conductance and current in unintentionally disordered double barrier GaAs- 
ALjjGai-zAs heterostructures. The source of disorder comes from interface 
roughness at the heteroj unctions (lateral disorder) as well as spatial inho- 
mogeneities of the Al mole fraction in the barriers (compositional disorder). 
Both lateral and compositional disorder break translational symmetry along 
the lateral direction and therefore electrons can be scattered off the growth 
direction. The model correctly describe channel mixing due to these elas- 
tic scattering events. In particular, for realistic degree of disorder, we have 
found that the effects of compositional disorder on transport properties are 
negligible as compared to the effects due to lateral disorder. 
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I. INTRODUCTION 



Resonant tunneling (RT) through double-barrier structures (DBS) make these systems 
very promising candidates for a new generation of ultra-high speed electronic devices. For 
instance, a GaAs-Gai_.jAL.As DBS operating at THz frequencies has been reported in the 
literature.0 The basic reason for RT to arise in DBS is a quantum phenomenon whose 
fundamental characteristics are by now well understood: There exists a dramatic increase 
of the electron transmission whenever the energy of the incident electron is close to one of 
the unoccupied quasi-bound-states inside the well.! In practice, a bias voltage is applied to 
shift the energy of this quasi-bound-state of nonzero width so that its center matches the 
Fermi level. Consequently, the j — V characteristics present negative differential resistance 
when the quasi-bound-state lies well below the Fermi level. 

In actual samples, however, the situation is much more complex than this simple pic- 
ture. Scattering by phonons, electrons or defects reduces the required quantum coherence 
and, in fact, deviations from the above simple description are observed. These scatter- 
ing mechanisms explain the occurrence of side-band resonant peaks due to the interaction 
with phononsiS and photons,iJi hysteresis in the j — V characteristics due to many body 
effects,^ and the decrease of electron mobility due to rough surfaces even in good-quality 
heterostructures.0 While one-dimensional models successfully describe scattering by elec- 
trons and phonons, they obviously cannot account for all interface roughness effects. So 
far, there are analytical results concerning the propagation of wave packets in a randomly 
layered medium when the potential is a random function of only one coordinate^ but for 
a small number of layers, as in DBS, in-plane disorder becomes important and one expects 
such approaches to fail. Realistic models of in-plane disorder, including power-like or ex- 
ponential spatial correlations observed by transmission electron microscopy, usually lead 
to intractable analytical models; hence the importance of numerically solvable models to 
bridge this gap. An important contribution was already provided by Henrickson et ai, who 
applied the tight-binding Green function method to interface roughness in DBS.0 In this 
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paper transmission was studied for a DBS with a rather simple model of disorder (periodic 
roughness with random relative phases at the interfaces). A level splitting was found for 
narrow well DBS and a dependence on the roughness lateral size was observed. 

In this work we present a simple two-dimensional model that allowed us to carry out an 
extensive study of electron scattering by interface roughness and how it manifests trough 
vertical transport properties of the DBS. In particular, we aim to elucidate the relationship 
between macroscopic properties of the DBS (e.g. conductance) with microscopic parame- 
ters of the model (e.g. correlation length of the surface disorder). The paper is organized 
as follows: The body of the paper is Sec. [0], where we present our model based, on the 
effective-mass approximation, describing a disordered sample (e.g. an imperfect DBS) con- 
nected to two perfectly ohmic leads. The Ben Daniel-Duke equation is discretized, boundary 
conditions are discussed and scattering solutions are found by means of the transfer-matrix 
method for any arbitrary heterostructure made of wide gap semiconductors. Afterwards, ex- 
pressions for two-terminal conductance and current for unintentionally disordered DBS are 
provided. The model is worked out in a two-dimensional space for computational limitations, 
although it will be clear that generalization to three dimensions is rather straightforward. 
Section |TTT] describes the various models of disorder we have used to mimic structural data 



obtained by transmission electron microscopy and X-ray scattering. In Sec. [TV] we present 
the numerical results and, in particular, the relationship between macroscopic properties of 
the DBS with microscopic parameters of the model. Finally, in Sec. [V] we discuss our results 
and the main conclusions of the work. 



II. MODEL 

A. Discrete equation 

We consider the single-electron two-dimensional Schrodinger equation within the frame- 
work of the one-band effective-mass approximation. Close to the V valley, this approach 
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leads to the Ben Daniel-Duke equation for the envelope function 

h 2 



+ — I + U(y,z) 



(1) 



2m* \dy 2 dz 2 

where z denotes the coordinate along the growth direction (see Fig. [I]). Notice that we have 
taken a constant effective mass m* at the T valley although this is not a serious limitation 
as the model can be easily generalized to include a position-dependent effective mass. We 
then consider a mesh with lattice spacings a y and a z in the y and z directions, respectively. 



Defining t y 



H /(2m*c) and t z = —% /(2m*a z ), we arrive at the following discretized 



version of Eq. ([!]) 



) + {U n ,m - 2t z - 2ty)lp n7m = Eljj 

n,i 



(2) 



The potential term U n>m in Eq. (|J) is given by the conduction-band edge energy at the 
point (na y ,ma z ) which, in turn, depends on the Al mole fraction in the vicinity of that 
position. Consequently, both kinds of disorder (lateral and compositional) are taken into 
account through this diagonal term (diagonal disorder). 

B. Transfer matrix formalism 



In order to solve the tight binding-like equation (§) we use the transfer matrix method 
based on the solutions calculated for each slide of the system along the z direction. To this 
end, we define 



n,2 



n = 0,1,...,N + 1, 



(3a) 



and 
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Here M and N + 1 are the number of mesh divisions in the y and z directions, respectively. 
Note that periodic boundary conditions have been taken into account on each slide. With 
these definitions, Eq. (0) can be cast in a more compact form 

\ / ^ \ 

(4) 
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t-\EI-Mn) ~1 

X O 



where X and O are the M x M identity and null matrices respectively. This expression 
allows us to relate by iteration the envelope function amplitudes at non-consecutive slides. 
In particular, we can obtain the expression of the envelope function amplitudes in the left 
contact as a function of the amplitudes in the right one 
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(5b) 



is the transfer matrix for the heterostructure. 



C. Scattering solutions 

We are interested in the scattering of electrons incident from the left contact. The 
envelope functions within the contacts will be determined by the boundary conditions. These 
boundary conditions are open in the z direction, and periodic on each slide, that is in the y 
direction. The former imply plane wave solutions in the z axis, and the latter yield an energy 
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discretization on y. As a consequence, this discretization results in a number of transverse 
channels equal to the number of points in the transverse mesh direction. Considering by the 
moment U n>m = at the contacts and setting ip n ^ = ip ny M+i, a particular solution of Eq. (0) 
is given by 

^n,m = ex P f * m J ex P ( ik i a z n ) > / = 1, 2, . . . , M. (6) 

Here TV; is a normalization constant that is needed for the different propagating modes to 
carry the same current all of them. This normalization constant is given by 

Ml = \ sin 2 [ H^j + -L sin 2 (kia z ) (7) 
o? y \ M J a? z 

Note that we consider a dependence on the index I of the momentum in the z direction. 
This is a consequence of the implicit assumption that all the scattering processes in region II 
are elastic, so energy is conserved. Thus, for a given energy 



ki = — cos < — 



( 2vrZ 
E - 2t v cos — - 1 



1 . (8) 



Finally we can write a general solution of the scattering problem for an electron of energy 
E impinging on the scattering region on the I transverse mode. On the left contact the 
solution will be of the form 

^n,m = 7^== ex P h -jf m I ex P {ikia z n) 

M 1 / 27T7 \ 

+ ^ f; j — exp ( % ———m J exp (—ikja z n) , (m, n) G I. (9a) 



3=1 

On the right contact we have 

M ~ 1 / 27T7 

Vim = XI *y "Ttt exp ( * ~TF m ) exp ( ik J a z n ) > ( m > n ) G IIL ( 9b ) 



IM 3 * V M 



3=1 

The matrices f and t appearing in these two solutions are the reflection and transmission 
matrices and they are responsible for the channel mixing due to scattering events. Thus, r^ 
represents the probability amplitude for an electron impinging in channel % to be reflected 
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into the chanel j. Similarly ty represents the probability amplitude for an electron impinging 
in chanel i to be transmited through the scattering region II into the channel j. Note that 
the solution within region II is unknown. Actually, we are not interested in this solution 
since all the physics of the scattering problem is contained in the mixing matrices t and f. 
In the following, our main goal will be to relate the elements of t and f to those of in 
Eq. ©. 

We start by re-writting formally Eq. (fj) as Hip = Eip. Now we perform the following 
transformation upon the envelope functions 



i/> = t~ V- 



(10) 



It is easy to see that the Hamiltonian is invariant under such a transformation, that is, TC = 
t^Tit = TC. This is clear since t and TC formally commute as they act upon different vector 
spaces. The fact that TC is invariant implies that T^' is not affected by transformation (|10"D 
so it remains unchanged. In addition, equations ( |9a|) and (|9T|) are transformed into 
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respectively, where we have defined bij = J2sLit ls r S j and aij = t 2 - . For each channel I we 
can write the transform of Eq. ([5a]) as 
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Now, introducing the definitions 
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we arrive at the following transform of Eq. ( |12D in terms of the elements of the mixing 
matrices 

(14) 
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and, finally, we have 



( tf\ ( j\o jgo\ 
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(15) 



/ 



that give us the matrix elements a~ij and by with j = 1,2, ... ,M for each channel / in terms of 
the product of two known matrices and the transformed solutions (|llb| ) at the right contact 
for the / channel. It is worth noticing that all the scattering information is contained in the 
transfer matrix T^ N \ Thus, this quantity turns out to be the fundamental object in the 
resolution of the scattering problem. 



D. Conductance and current 

Once a and 6 have been calculated, obtaining t and r is an easy task. However, the 
matrices calculated this way are the response matrices that contains some non-physical in- 
formation. This non-physical information comes from the fact that we have been considering 
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all the mathematical solutions of the scattering problem, including those diverging at infin- 
ity. These solutions are those for which momentum in Eq. @ is an imaginary number. In 
order to avoid the unphysical solutions we will cancel out all the elements on t for which the 
incoming and/or the outgoing transverse quantum number I satisfy the following condition 

E < 2ty [cos — - lj . (16) 

Once the physical mixing matrices are known, particularly the transmission matrix t, we can 
use them to calculate different physical quantities like the conductance or the electric cur- 
rent. From the Landauer-Biittiker formalismjl! the zero temperature two-leads multichannel 
conductance can be calculated using de Fisher-Lee formula0 

G =^-Tr(Vt). (17) 

In order to calculate the electric current due to an applied field F we have to modify slightly 
our equations and further approximations are needed. First of all, we consider perfect leads, 
so no voltage drop occurs within the contacts. Besides, we will not consider the contribution 
to the current of electrons impinging from the right contact. This is a good approximation 
provided eV > Ep, where V is the voltage drop along the region I and Ep is the electron 
Fermi level on the contacts. Finally, we will assume zero temperature so no carrier statistics 
will be taken into account. 

When an electric field applied along the growth direction is considered, the potential 
U(y,z) in Eq. (U) has to be replaced by U(y,z) + Up(z), where Up is constant in the 
contacts (Up(z) = in region I and Up(z) = —eV in region III) and behaves linearly in 
region II, namely Up(z) = —eFz. Here F is the applied electric field and V = FL is the 
voltage drop across the scattering region whose length is L. The electron momentum on the 
right contact now reads 



1 if 1 

qi = — cos i — 

a z { 2t-, 

and solution in region III changes to 



/ 2vr/ 
E + eV - 2L, cos — - 1 



M 



l\, (m,n)e III (18) 



M ~ 1 / 2nj 

i> l n, m = E h -777 exp r ~M m ) exp (*3? a * n ) • ( 19 ) 



Bearing in mind all the approximations considered so far, a new transfer matrix 
can be calculated for a each value of F. Inserting solution (|T9| ) in Eq. ( |T5| ) we can obtain 
the mixing matrices for a given applied field and then calculate using them the electric 
current induced by that field. From the expression for the probability current we can write 
the discretized version of the electronic current density contribution of chanel I along the 
applied field direction across the right lead (calculated on a slide n within the region III) as 

pf n 1 M 

A = i J £rT7 E [« M + i,m - < m (^ +1 ,J*] • (20) 
na y M m=i 

In order to calculate the total density current across the sample, we just need to sum up 
over all the allowed transverse states below the Fermi level at the left contact, that is 

Notice that the electron momentum in the z direction, k, varies continously as the energy 
for an incident electron can take values from to Ep. Here, k l F is the component of the 
electron momentum along the z direction for an incident electron with transverse quantum 
number I and an incident energy Ep. Using the transmission matrix elements we finally 
arrive at 



et z a z 1 M M r k F r l # l l l -i 

3z — ^ 2-nfia ~M ^ ^ Jo ( ^71,771 ) ^n+i.m — (^71+1,771) (21) 



et z a z 1 MM kl MM^^ l / 2 n(s-j] 



-m 



x exp [i(q s — qj)a z n] [exp (iq s a z ) — exp (—iqja z )] dk. (22) 

The results we have obtained so far provide an exact, although nonclosed, analyti- 
cal description of any two-dimensional heterostructure based on wide-gap semiconductors, 
whenever the Ben Daniel-Duke holds valid. Let us stress that the generalization to three- 
dimensional heterostructures is fairly straightforward. With these results at hand, we can 
compute the transport magnitudes we mentioned above. All expressions are very simple 
and suitable for an efficient numerical treatment for any specific case. We will now evaluate 
them for unintentionally disordered DBS to describe the relevant features of the model. 
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III. MODELS OF DISORDER 



Unintentional disorder appearing during growth in actual heterostructures depends crit- 
ically on the growth conditions. There exist several techniques, like scanning tunneling 
microscopylil~0 and X-ray scattering,@ which have been applied in recent years to quanti- 
tatively determine structural properties of multilayers and superlattices. Precise information 
about the nature and extend of defects at interfaces is now available. Following Mader et 
ai, disorder in a heterostructure can be classified into two categories, namely lateral and 
vertical.0 Lateral disorder occurs whenever one semiconductor protrudes into the other, 
forming chemically intermixed interfaces, steps and islands. As a consequence, the interface 
is not flat and translational symmetry in the plane perpendicular to the growth direction is 
broken. On the other side, vertical disorder is observed whenever layer thicknesses fluctuate 
around their nominal values. Compositional disorder due to different local Al concentration 
in each Gai_ x ALAs layer can be viewed as mixing of vertical and lateral disorder. 



A. Lateral disorder 

In order to treat lateral disorder we have considered the formation of islands on the 
interface between two consecutive layers having identical lateral sizes all of them and being 
consecutive one to each other. In our model islands have heights that are randomly dis- 
tributed (see Fig. 0). It is possible to express the rough profile of the interface between two 
consecutive layers defining the following height function 

%) = t) £ w n{Q(y - nC) + Q[(n + 1)C -y)- 1}. (23) 

n 

Here h(y) represents the deviation from the flat surface at position y, is the Heavyside 
theta function, ( is the island width, w n is a random variable associated to the nth island 
that controls the fluctuation around the mean value, and r\ is the largest deviation - 
in absolute value — assuming that the w„'s are uniformly distributed between —1 and 1. 
Hereafter rj will be referred to as degree of lateral disorder. In the following, two models 
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will be introduced. First of all, uncorrelated disorder will be considered, for which the 
random variables w n take values from —1 to 1 uniformly, satisfying the following correlator 
(w n w m ) = S nm /3. However, experimental data and theoretical models of molecular-beam- 
epitaxy growth processes^ indicate that the heigths of the surface at distant points are 
strongly correlated. To mimic this situation, we will consider a second model of lateral 
disorder where the w n 's are normally distributed with zero mean satisfying the exponential 
correlator (w n w m ) = exp(— \n — m|/£)/3, where £ is the correlation length of the disorder. 

B. Compositional disorder 

Compositional disorder in GaAs-Al x Gai_a;As heterostructures is due to the lack of spatial 
uniformity of the Al mole fraction during the heterostructure growth process. We will simu- 
late compositional disorder by defining a spatial two-dimensional mesh over the Al x Gai_ x As 
barriers, as shown in Fig. |^. We will associate a value for the Al mole fraction to each 
region of the mesh. This value will fluctuate randomly around the nominal value of the 
aluminum fraction on the barrier x. Thus we can write 

Xij = x + flwij (24) 

where Q is the maximum fluctuation of the Al mole fraction, and Wij is a random variable 
related to the point of the mesh, taking values uniformly between —1 and 1. Note that 
the overall averaged mole fraction is equal to x. 

IV. RESULTS 

We have performed several numerical calculations in order to study the effect of both 
lateral and compositional disorder over the transport properties of DBS made of GaAs- 
Al x Gai_ x As heterostructures. 
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A. Uncorrelated lateral disorder 



We start by considering the effect of the interface roughness, characterized by the degree 
of lateral disorder r\ given in fl23|), on the conductance. Figure |] shows the conductance 
calculated for different values of rj. Here a y = 10 nm, a z = 0.3 nm, ( = 20 nm, M = 50, 
and N = 38. The barrier widths are 2.1 nm for both the emitter and the collector, their 
heights are also the same, 0.3 eV, and the well width is 4.8 nm. The curves correspond to 
an ensemble average of the conductance curves over 100 different realizations of the lateral 
disorder. Three different values of r\ were studied, namely rj = (perfect DBS), rj = 0.3 nm 
(largest fluctuation of the order of one monolayer) and rj = 0.6 nm (largest fluctuation of 
the order of two monolayers). As a main result, it can be seen in Fig. |] that increasing the 
degree of lateral disorder, rj, results in a decrease of the conductance at the resonant energy. 

Notice that the resonant conductance peak slightly widens due to the fluctuations of its 
energy for each realization of the disorder. Besides, an additional effect can be seen, that is, 
as the degree of lateral disorder rj increases the conductance peak shifts to smaller energies. 
This effect will be important in order to explain the j — V characteristics later. Regarding 
the effect of the size of the islands, (, we have observed that it can be neglected unless this 
size is of the order of the electron wavelength, that is, for £ ^> A e the conductance does not 
depend on (. As expected, when ( ~ A e the electron starts to see the disorder and then 
the conductance increases as £ decreases, as shown in Fig. ||. For energies about 0.1 eV this 
transition takes place at sizes £ ~ lOnmS 

Current as a function of the applied bias is depicted in Fig. [| For this calculation we 
have chosen a y = 10 nm, a z = 0.1 nm, ( = 10 nm, M = 20 and N = 77. The barrier 
widths are 2.1 nm for both the emitter and the collector and their heights are also the same, 
0.25 eV. The well width is 2.9 nm. Curves show in Fig. ^ correspond to an average over 
50 different realizations of the lateral disorder. The Fermi level was fixed at 20 meV, the 
maximum applied bias was 0.5 eV, and r\ = 0.3 nm. Surprisingly, current in disordered DBS 
is larger than that measured in ordered DBS. This counterintuitive result can be explained 
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recalling that conductance peak shifted to lower energies, as mentioned above. The lowering 
of the conductance peak can be understood assuming that surface roughness makes the 
effective width of the quantum well larger than its nominal value. This would imply a lower 
transmission resonance so the current peak shifts toward lower bias. But an effective wider 
well implies a higher current too, so in some statistical sense the current for the disordered 
DBS will be higher. 

B. Correlated lateral disorder 

Up to now we have been considering only lateral uncorrelated roughness at the hetero- 
junctions. We have also calculated the conductance for a DBS in which interface roughness 
is an exponentially correlated random variable, as described previously. Figure [7| shows the 
conductance at a given energy about the conductance resonance for the ordered case (where 
fluctuations are lower), E c = 0.1025 eV, when the correlation length varies several orders of 
magnitude. Physical parameters of the DBS are the same as in Fig. |j. The size of the islands 
is 20 nm and the degree of lateral disorder is 77 = 0.3 nm. An average over 50 realizations 
of the disorder has been performed. It can be observed that the conductance increases as 
the correlation length increases. Two asymptotic regimes are clearly observed for small and 
large correlation lengths, respectively. The former limit, £ —>■ 0, can be understood as the 
situation previously studied where lateral disorder is uncorrelated. From a closed inspection 
of Fig. [7] we come to the remarkable result that this uncorrelated limit correctly describes 
transport properties whenever the correlation length do not exceed 10 nm. Then we are led 
to an important conclusion, namely those models of electron transport in disordered DBS 
based on the assumption that disorder is uncorrelated would yield right values of the con- 
ductance even if the correlation length is not vanishing but fairly large. There is no need to 
mention that uncorrelated disordered models are much more convenient for analytical work 
than correlated ones. The size of the plateau for which correlations in the disorder do not 
play an essential role is governed by the size of the islands, (, as might be expected. This 
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is due to the fact that only correlations of the lateral disorder whose correlation length is 
greater than the size of the islands affect electron motion. For these sizes, the conductance 
increases as the correlation length becomes larger. This trend is observed up to sizes of the 
order of 100 nm, namely the system size in the y direction, and then the second regime ap- 
pears. In this regime the limiting value of the conductance can be calculated as an average 
over different ordered DBS, each one having barrier widths distributed normally. 

C. Compositional disorder 

Having discussed the effects of interface roughness, let us now present our results re- 
garding compositional disorder. To clarify the effects of local fluctuations of the Al mole 
fraction, we assume that the interfaces are perfectly flat and that the width of the different 
semiconductor layers are exactly their nominal values, i.e., we neglect fluctuations of their 
widths. We show typical results in Fig. |8|, where physical parameters are the same as in 
Fig. U The espatial extend of the region where the Al mole fraction can be assumed locally 
constant was taken to be a = 10 nm and b = lnm (see Fig. |^). It can be observed that, for 
fluctuations as large as Q = 0.14, that is, height fluctuations of about 0.1 eV, when aver- 
aging over 100 realizations the conductance is almost unchanged (only 3% variation). This 
behavior is a consequence of the small electron probability amplitudes within the barriers, 
where compositional disorder is considered to appear. In other words, even moderately high 
spatial inhomogeneities in the barriers scarcely affect the electron envelope function. Thus 
we can conclude that the effects of compositional disorder can be disregarded as they are 
much smaller than those of lateral disorder. 

V. CONCLUSIONS 

In this paper we have presented a method to study electron transport in unintention- 
ally disordered DBS when translational symmetry along the normal plane to the growth 
direction is broken by structural imperfections. Two different kinds of disorder leading to 

15 



electron scattering has been presented, namely lateral disorder and compositional disorder. 
Lateral disorder appears as a consequence of the roughness present at the interfaces between 
GaAs and ALjGai-^As. Compositional disorder is due to the inhomogeneity of the Al mole 
fraction occurring during the growth process of ALjGai-aAs epilayers. We have shown that 
the main effect of the lateral disorder is to decrease the DBS conductance. When correla- 
tions are introduced, conductance starts improving, the larger the correlation length, the 
higher the conductance. However, current shows a surprising behavior, as it is higher in 
the disordered case than in the ordered one. We attributed it to an energy shift of the 
transmission resonance of the DBS. Another important conclusion regards the validity of 
uncorrelated disordered models; they provide right values of the conductance whenever the 
actual correlation length does not exceed the system size, thus allowing a simpler theoret- 
ical description. Finally we have shown that the effect of compositional disorder can be 
disregarded as it is much smaller than the effect of lateral disorder. 
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FIGURES 

FIG. 1. Schematic view of the sample. Regions I and III are the electrical leads of the sample 
(contacts) and electrons undergo scattering processes only at region II. 

FIG. 2. Steps and islands modeling lateral disorder (roughness) at the interface between two 
epilayers. 

FIG. 3. Schematic view of the Al fraction mesh in the barrier showing the distribution of Al 
concentration at different sites. 

FIG. 4. Conductance in units of 2e 2 /h of an ordered DBS compared with that of a disordered 
DBS for two different values of the degree of lateral disorder r\. The disordered results were obtained 
by averaging over 100 realizations of the disorder. The inset shows the shift of the first resonant 
peak towards lower energy on increasing the degree of disorder. 

FIG. 5. Current in arbitrary units as a function of the applied bias along the z direction for an 
ordered DBS and a disordered one. Curve was calculated by averaging over 50 realizations of the 
disorder. 

FIG. 6. Conductance in units of 2e 2 /h as a function of the size of the islands for ( ~ A e (dotted 
line) and £ < A e (dashed line). Curves were calculated by averaging over 50 realizations of the 
disorder. 

FIG. 7. Conductance in units of 2e 2 /h of a disordered DBS where disorder is exponentially 
correlated, as a function of the correlation length, for a fixed energy E c = 102.5 meV. Each point 
of the curve was calculated by averaging over 50 realizations of the correlated disorder with fixed 
correlation length £. 

FIG. 8. Conductance in units of 2e 2 /h for an ordered DBS compared with that of a DBS with 
only compositional disorder for two different values of Q. The inset shows enlarged view of the 
first resonant peak. 
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